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On  pane  14,  the  integrand  of  the  integral  / should  be  multiplied  by 

the  product  of  tho  scale  factors  h^(£o,t])  ant'  S /'0,T])  • ^Rua^ion 
now  reads: 

h j = i S fnj2 

m'"  -nji  “S 

n ,2 

Tho  integral  / S (r|)dr,  in  Eqs.  3.2  and  3.3  is  now  replaced  by  the 

-Si 

integral : 

(if  (<£")* jfot  A.W- 

nji  fiz  Cijs 

It  is  seen  that  tho  product  of  the  two  integrals  < and  j can 

^il  ^jl 

now  be  written  oxplicitly  as  follows,  after  substitution  of  Eq.  6.4: 

“Si 


luiwwi  lor  pobUo  rttooNt 
OymbofloB  Unlimited 


RADIATION  IMPEDANCE  FOR  THE  ELEMENTS  OF  A 


V 


CYLINDRICAL  TRANSDUCER  IDEALIZED  AS  AN  efeLATE  SPHEROID. 


A/jechnical  jAepmt, 


Preoarod 


f or 


The  Department  of  the  Navy 
Office  of  NavaiJiflSflaLC 
Under  Contract/  Nonr-2739 


ice  of  nlaVfi 
ContracyT 

(0- 


by 


l.fli 


Richard  V, /Waterhouse 


CA"BR  IDGE  ACQUST  ,CAL  ASSOCIATES,  INC. 
1278  Massachusetts  Avenue 
Cambridge  313,  Massachusetts 


0 7 7~  ?6~0 


■ 


TABLE  OF  CONTE  MI'S 

Page 

1 Introduction  ..o..u...< 1 

2 Comparison  of  Transducer  Element  and  Surface  Element 

of  Spheroid  2 

3 Choice  of  Spheroid  to  Represent  Transducer  4 

4 Oblate  Spheroidal  Coordinates  9 

5 Solution  of  the  Wave  Equation  in  Spheroidal  Coordinates  ...  11 

6 Boundary  Cond it  ion 11 

7 Near-Field  Pressure  14 

8 /Mutual  Impedance  Between  Two  Elements 14 

9 Degeneracy  of  Z 15 

^ J 

10  Computation  Required 17 

List  of  Symbols 20 


ACCESSION  for 

HT|S  White  Section 

POC  Bull  Section 


I UNANNOUNCED 

HEZfeB r4 


BY — — 

mtIBIlTWIfMMUBlUn  COOES 


KL*  avail  ind/or  st-tu/wl 

t 

j 

1 


□ □ 


RADIATION  IMPEDAiJCE  F(Jfs  THE  ELEMENTS  OF  A 
CYLINDRICAL  TRANSDUCER  IDEALIZED  AS  AN  OBLATE  SPHEROID 


1 . Introduction 

FThe  transducer  in  question  has  the  form  of  a ripht  circular 
cylinder;  the  radiating  elements  of  the  transducer  correspond  to  the  curved 
surface  of  the  cylinder,  and  the  plane  surfaces  of  the  cylinder  correspond 
to  the  top  and  bottom  of  the  transducer,  which  are  covered  by  rigid  plates 
and  are  assumed  to  have  no  motion.  The  transducer  elements  have  square, 
plane  faces  and  form  a mosaic;  the  total  radiating  surface  comprises  several 
hundred  elements.  Each  element  can  be  driven  separately,  and  the  problem 
is  to  calculate  the  mutual  radiation  impedance  between  any  two  elements, 
situated  at  arbitrary  positions.  From  this  information  the  impedance  when 
more  than  one  element  is  driven  can  be  found  easily  by  superpos i t ion.  'f — \ 
The  problem  is  one  of  a large  group  of  problems  which  consist  of 
Calculation  the  mutual  impedance  of  two  sources  placed  in  various  environ- 
ments. For  example,  the  sources  may  be  (a)  in  a free  field, ^ (b)  in  the 
half-space  above  a rigid  reflecting  plane,  (c)  on  the  surface  of  such  a 

plane,  ’ (d)  on  the  surface  of  a rinid  sphere,'*  (e)  on  the  surface  ot  a 

6 7 

rigid  infinite  cylinder.  * 

The  sources  may  be  of  monopole,  dipole,  piston  or  other  typo,  and 
thus  the  number  of  problems  involved  is  quite  large;  the  references  cited 
give  some  results  for  those  problems  that  have  appeared  in  the  literature. 

The  present  problem  involves  calculating  the  near  sound  field  of 
the  cylindrical  transducer  when  any  one  element  is  drivon.  However,  there 

exist  no  exact  solutions  for  the  near  sound  fiold  of  a finite  cylinder,  and 

6  8 

the  approx  invite  solutions  * involve  mathemat ical  forms  whose  computation 
would  require  a very  large  effort. 

1 Rayleiph,  Collected  Papers  (Cambridpe  University  Press,  London,  1912), 
Vol  5,  p.  137 

2 S.  J0  Klapman,  J.  Acnust.  Soc.  Am.  1 1 . 2B9,  1940 

3 R.  L.  Pritchard,  Tech.  Memo.  No.  21,  Appendix  C,  Acoustics  Research 
Laboratory,  Harvard  University,  NR-01 4-903  (Jan.  1b,  1951) 

4 R.  V.  Waterhouse,  J.  Acoust.  Soc.  Am.  30,  4,  (1958) 

5 C.  H.  Sherman,  J.  Acoust.  Soc.  Am.  31 . 947,  (1959) 

6 D.  T.  Laird  and  H.  Cohen,  J.  Acoust.  Soc.  Am.  24.  46,  1952 

7 D H.  Robey,  J.  Acoust.  Soc.  Am.  22,  706,  (1955) 

8 W.  E.  Williams,  Proc.  Camb.  Phil.  Soc.  Jj2 , Part  2,  322,  1956 
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It  is  therefore  considered  more  expedient  to  represent  the 
cylindrical  transducer  as  an  oblate  spheroid,  in  which  coordinate  system 
•an  exact  solution  for  the  near  field  is  available,  althou  h the  functions 
involved  are  not  well  tabulated.  This  representation  by  a spheroid  is  ex- 
pected to  give  more  accur  ite  results  than  the  representation  by  a sohore 
which  has  been  studied  by  Sherman  at  the  U.  S.  Navy  Underwater  Sound 
Laboratory  at  ew  London,  Connecticut. 

Thus  the  problem  to  be  solved  is  the  boundary  value  one  of  a 
quasi-roctangular  source  situated  on  the  surface  of  a rigid  oblate  spheroid; 
as  far  as  is  known  this  problem  has  not  been  treated  before.  The  source 
configuration  is  described  as  quasi-rectangular  as  it  is  the  portion  of 
surface  of  the  spheroid  bounded  by  the  orthogonal  coordinate  lines;  this 
surface  approx imatos,  but  is  not  identical  to,  a plane  rectangle  as  discussed 
further  in  Section  2. 

2.  Comparison  of  transducer, -Clement  and  Surface  ' lemon t of  Spheroid 


The  transducer  is  shown  schematically  in  Fig.  1.  Its  height  is 
3.9\  3nd  its  diameter  11. 4\;  there  are  IxG^AlS  mosaic  elements,  each  of 
which  has  a W2  square  radiating  face. 

A surface  element  on  the  spheroid,  Av|A0,  is  bounded  equator  ially 
by  the  parallel,  curved,  coordinate  1 ines  r(  where  h|ji“y)j2=  and 

vertically  by  the  nonparallel,  curved,  coordinate  lines  0.^,0. where 
0 i ^ “0 i 2 = The  surface  elernont  is  not  plane.  Such  an  element  will  bo 

used  to  represent  the  plane,  square  surface  of  a transducer  element,  see  Fig.  5. 

There  are  seven  rows  of  elements  comprising  the  transducer;  i 

however,  by  symmetry  only  the  four  elements  numbered  1 to  A in  Fig.  1 are 
needed  to  characterize  the  pressure  field.  The  pressure  field  over  the 
spheroid  will  be  calculated  for  each  of  tho  four  elements  driven  separately. 

There  are  three  main  differences  between  the  spheroid  elements 
At)A0  and  tho  transducer  elements  they  represent,  and  those  will  bo  discussed 
in  turn.  These  differences  are  greatest  for  the  top  row  of  elements,  con- 
taining element  1 in  Fig.  1. 

a.  The  spheroid  elements  exceed  the  transducer  elements  in 


surface  area,  the  largest  difference  being  in  the  top  ro w where  for J =0.60 
the  areas  are  in  the  ratio  of  approximately  1.8:1;  for j =0.30  the  ratio  is 
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approximately  1.4  to  1.  Thir.  area  difference  will  affect  the  results  in 
which  the  pressure  is  integrated  over  the  surface  lenient,  as  in  8.1  „ It 
would  bo  possible  to  compens >te  for  this  area  difference  by  taking  the 
spheroid  surface  for  position  1 to  ho  less  than  the  projected  area  of  the 
transducer  element,  i.e.,  taking  the  annular  coord inate»'| ^ correspond ing 
to  the  transducer  edge  to  be  slightly  closer  to  the  equator  of  Ihn  spheroid 
than  before.  However,  this  change  would  also  make  the  element  less  square, 
and  more  elongated,  and  would  affect  the  boundary  condition  (6.1).  Thus  it 
is  probably  better  not  to  make  this  correction. 

b.  The  nonparallel  coordinate  lines  0 - ^ j 2 r®Prf>s^n^  tho  vertical, 

parallel  boundaries  of  the  transducer  element,  and  the  nonparallel  ism  is 
greatest  for  the  top  row  of  elements. 

c.  The  coordinates  1 ines  ^ representing  the  longitudinal 

edges  of  the  transducer  element  are  unequal  in  length,  the  difference  is 
greatest  for  the  top  row  of  elements,  where  it  is  approximately  13  °/o  for 
> 0.60. 

The  above  differences  between  the  element  of  spheroid  surface  and 
the  transducer  element  it  represents  will  cause  the  results  to  deviate 
slightly  from  those  applyinp  to  a cylindrical  transducer.  However,  this 
is  the  price  paid  for  the  relative  convenience  of  using  this  coordinate 
system  for  which  exact  solutions  to  the  wave  equation  can  be  written  down. 

3.  ChP-ico  p.f.Jjphu,r.ai.£Ll,p  ^argsant,  Transducer 

In  choosing  the  spheroid  that  will  best  represent  the  circular 
cylindor,  a compromise  has  to  be  made.  Mgs.  1 and  2 show  the  ellipses  and 
rectangle  which  are  the  cross  sections  of  the  sphoroicfeand  cylinder; 
rotition  about  the  z axis  gives  the  3-dimensional  figures.  The  equatorial 
portion  of  the  spheroid  corresponding  to  the  vibrating  surface  of  the 
transducer  is  that  contained  within  the  two  long  sides  of  the  rectangle. 

If  the  spheroid  chosen  has  a polar  diameter  nearly  equal  to  the  height  of 
the  transducer,  the  portion  of  the  sphoroid  corresponding  to  the  radiating 
surface  of  the  transducer  is  too  large  and  vory  curved,  and  the  normal 
motion  of  the  sphoroid  element  will  not  represent  closely  the  normal  motion 
of  the  cylindor  element. 
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However,  if  the  polar  diameter  of  the  spheroid  is  made  large,  so 
that  the  part  corresponding  to  the  radiating  surface  is  more  nearly  parallel 
to  the  latter,  the  discrepancy  between  polar  diameter  and  cylinder  height 
becomes  larger;  also  in  this  case  the  spheroid  approaches  a sphere,  for  which 
the  corresponding  results  are  already  known.  Thus  a compromise  must  be  made 
between  the  fit  at  the  vibrating  surface  and  the  fit  at  the  onds  of  the 
transducer. 

A factor  bearing  or  the  degree  of  nonoaral lei  ism  tolerable  between 

the  vibrating  surface  of  the  cylinder  and  the  corresponding  surface  of  the 

spheroid,  is  the  directional  pattern  of  the  transducer  element.  In  the 

oresent  case  t h* * element  is  X/2  square,  and  +hus  its  far  field  radiation  is 

9 

not  highly  directional,  as  shown  in  Fig.  3.  This  implies  that  the  results 
will  not  be  sensitive  to  this  nonpiral lei  ism,  as  they  would  bo  if  the 
elements  wore  highly  directional. 

In  the  lioht  of  all  these  factors,  the  spheroid  whose  proportions 
are  specified  by  the  value  jQ=f .6  is  considered  a satisfactory  compromise, 
and  this  value  is  used  in  J ho  followino  discussion.  However,  another  value 
of  3o  could  be  chosen  without  causing  any  significant  changes  in  the  mathe- 
matical aooroach  in  this  resort.  The  spheroid  corresDond inp  to  5 =0„6  is 

^ o 

shown  in  Fig.  1;  its  major/minor  axes  are  in  the  ratio  1.9  to  1.  It  thus 
differs  aopreciably  from  a sphere.  The  minor  axis  of  the  spheroid  = 1.61  x the 
height  of  the  transducer,  and  the  normals  to  the  spheroid  surface  at  the 
points  corresponding  to  the  edges  of  the  transducer  are  inclined  at  - 57°  to 
the  x axis,  see  Fig.  4. 

The  length  of  the  rectanglo  was  taken  as  0.939  of  the  long  axis 
of  the  sph-roidj  this  was  considered  to  give  a better  match  at  the  vibrating 
surface  than  taking  it  the  same  length  as  the  long  axis.  This  value  tends 
to  minimize  the  average  difference  between  the  separation  of  two  elements 
on  the  transducer  and  on  the  spheroid.  The  values  of  the  relevant  parameters 
for  spheroids  of  various  proportions  are  given  in  Table  1,  and  Fig.  2 shows 
sections  of  sphoroids  corresponding  to  5o=0»4  and  0.8.  If  the  problem  is 
of  sufficient  interost  to  justify  the  computing  time  needed,  results  could 
be  computed  for  more  than  ono  spheroid. 


• » 


$ H.  Olson,  Elements  of  Acoustical  Engineering  (Van  Nostrand  Co0,  N.  V 
2nd  edition)  p. 32 


Diroctioru.il  pattern  of  a plane,  square  piston 
vibrating  in  an  infinite  rigid  baffle.  Graph 
shows  relative  pressure  as  a function  of  polar 
an^lo  in  a plane  parallel  to  tho  direction  of 
motion  of  the  piston  and  parallel  to  one  of  its 
sides.  See  Ref.  9. 
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TABLE  1 


T 

* 

Frequency 
parameter  ; 
c 

Polar  d jams  ter. 

transducer  height 

i 

Inclination  to 

x axis  of  normal  at 
edge  of  transducer 

~~.4  ’ 

35.6 

1.16 

78° 

.5 

34.2 

1.39 

66° 

.6 

32.8 

1.61 

57° 

.7 

31.3 

1.79 

50° 

.8 

29.9 

1.95 

44 J 

.9 

28.4 

2.08 

39° 

1.0  | 

27.1 

2.21 

1 35° 

Paramo tors  of  various  oblate  spheroids  which 
could  bo  usod  to  represent  the  cylindrical  transducer 


4 . Oblate  Spheroidal  Coordinator 


These  coord  incites  are  obtained  from  elliptical  coordinates  by 
rotating  them  about  the  minor  axis  of  the  ellipse.  They  are  obtained  from 
rectangular  coordinates  by  the  transf armat ion 


x=|[(1-»,2)(52+1)]^cos0 


y = 4 [(1-rj2)Cj2+1  )]a  Sin  0 


2 

= 4 -an 


• 2 


^5 


(4.1) 

(4.2) 

(4.3) 


• ..„,r  _ \ ..  > <r_ 


w 


The  coordinate  system  is  represented  in  Fig.  5;  the  symbols  are 


defined  in  the  list  of  symbols,  and  the  notation  follows  that  of  Flammor. 
The  surface  *=  constant  is  the  surface  of  an  oblate  spheroid;  the  surface 
constant  is  a rectangular  hyperboloid  of  two  sheets;  the  surface 
0 = constant  is  a plane  containing  the  z axis,  see  Fig.  5. 


10 


10  C.  Flammor,  "Spheroidal  Wave  Functions,"  Stanford  Un i vers i ty  Press, 
Stanford,  Cal.  (1957; 
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5.  Solution  of.  W^ye  Equation  in  Spho.ro id 3 1 Coordinates 

Oblate  spheroidal  coordinates  are  one  of  the  eleven  coordinate 
systems  in  which  the  wave  equation  is  separable.  When  sinusoidal  time- 
dependence  is  assumed,  the  wave  equation  reduces  to  the  Helmholtz  equation, 
whose  solution  in  oblate  spheroidal  coordinates  for  diverging  waves  is 


f &V0>  - Y. 

mf  n 


(5.1) 


The  symbols  are  defined  in  l.ho  list  of  symbols;  Hv*  not  it  ion 
follows  that  of  Flammer . S (-ic,>i)  is  an  annular  function;  and  R^^(-ic,  if) 
is  a radial  function,  -‘e  can  neglect  the  annular  functions  of  the 
2 kind,  Smn  , which  tend  to  infinity  at  the  poles, >)  = i 1,  where  we 
require  to  be  finite.  All  the  angular  functions  of  r|  used  here  will  be 
of  the  first  kind,  and  will  be  written  5 (n),  without  superscript.  In  this 

problem  tho  frequency  parameter  c is  fixed,  so  we  will  drop  it  from  tho 
exhibited  variables  ('Vj,0,^)  on  which  the  radial  and  angular  functions  depend. 

Tho  coefficients  amnj  in  5.1  must  now  be  calculated  from  the  main 
boundary  condition. 


6.  Boundary  Condition 


The  main  boundary  condition  is  tho  Neumann  one  in  which  the 
normal  derivative  of  the  dependent  variable  in  the  wave  equation  is 
soecified  over  an  interior  surface;  i.e.,  the  particle  velocity  is  specified 
all  over  the  spheroid. 

Thus  on  the  spheroid  surface  the  given  or  prescribed  normal 

velocity  is 


vg/>^ 

th 


= 1 


f fail 

1 1C0.  v<  K< 

ewhore 


(6.1) 


cl s< 


This  corresponds  to  tho  i tr  risducer  element,  at  tho  angular  position 
V-Ua'-W.  being  driven  with  velocity  u while  tho  rest  of  the 
surface  remains  motionless.  The  element  has  angular  dimensions  20^  in 
longitude,  *n  lJ^*tude,  v”here  i=1  to  A specifies  which  of  tho 
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four  transducer  elements  i3  excited. 

without  loss  of  generality  the  elomont  can  be  taken  as  centered 
on  the  coordinate  line  6^0  on  the  spheroid.  Then  since  the  boundary  con- 
dition requires  non-zero  velocity  at  {5=0,  the  solutions  in  5.1  containing 
sin  m0  are  eliminated. 

From  5.1  the  velocity  normal  to  the  soheroid  is 


vr  O|,0)  = 


1 


& 


'5 


amni  Rmn),(^}  Smn(>l)c0s  m 0 


(6.2) 

(6.3) 


m,  n 


In  6.2,  hv(-j)  is  the  scale  factor  applying  to  tho  coordinate  system  and  is 
given  by 

w*«(W 

where  d is  the  interfocal  distance  of  the  spheroid. 

Vq  need  an  expansion  for  the  given  velocity  (6.1)  in  a series  to 
match  6.3.  However,  a complication  arises  because  the  scale  factor  h^) 
in  6.3  is  a function  ofy^ 
of  sohoroidal  harmonics 

v 


(6.4) 


25 


If  the  given  velocity  is  expanded  in  a series 

M>  = ^_imni  Smn(V  cos  m 0 (6*4a) 


m,n 


the  coeff ic ients . cannot  be  matched  with  the  a • in  6.3,  since  the 
mn  i mn l 

litter  contains  hjtv-  An  expansion  of  the  given  velocity  in  terms  of  the 
soheroidal  harmonics  6.4a  divided  by  hs(’))  is  unfruitful,  as  the  resulting 
functions  do  not  form  an  orthogonal  set. 

The  difficulty  can  be  circumvented  by  expanding  tho  product 


h ^ V(>0)  = Ijmni  Smn(V  c05  m 0 

m,n 


(6.5) 


and  dividing  out  tho  scale  factor  later,  (see  6. 9)u  The  coefficients  4ni 
are  then  obtained  as  usual  by  multiplying  (6.5)  by  ^ (»j)  cos  m 0 and 

integrating.  Vie  define  an  integral 


jJ 


X. 


th  th 

This  gives  the  mutual  impedance  between  the  i and  j element* 


of  the  transducer.  The  calculation  of  these  mutual  impedances  is  discussed 


in  Section  10.  The  number  of  different  values  of  Z . needed  to  specify  this 


quantity  for  every  element  of  the  transducer  is  544,  as  discussed  in  Section  9. 
The  total  radiation  impedance  imposed  on  the  element  when  an 


arbitrary  number  N of  other  elements  are  being  driven  is  found  by  superposition 
to  be  N 

zj-^Zzu“«  <!M) 

J / = 1 


where  the  i source  is  driven  with  velocity  Uj 


9.  Degeneracy  cf  l , 


The  transducer  consists  of  7 rows  and  68  columns  of  elements, 
comprising  476  elements  in  all.  The  ouestion  arises  as  to  how  many  different 


values  aro  needed  to  specify  the  mutual  impedance  Z j j between  every  pair  of 


elements  i and  j. 

First  we  note  that  by  reciprocity 


Z.  . = Z .. 

•J  J* 


(9.1) 


This  reduces  the  number  of  different  values  of  l j .,  but  even  so  if  the  trans- 

J 

ducer  oossessed  no  spatial  symmetry  of  any  kind,  the  number  of  different 


values  would  be  476  (1  + 476/2)  = 113,764,  a large  number. 

In  fact  the  transducer  has  rotational  symmetry  about  its  axis  and 


Table  2. 


source  o 

laments  and 

any  column 

zu 

• 

• 

TABLE  2 

• • 

Csi 

5T" 

Z22 

• 

• • 

Z13 

Z23 

Z33 

• • 

ZU 

Z24 

Z34 

Z44  ' 

Zl5 

Z25 

Z35 

• • 

Z16 

Z26 

• 

• • 

Z17 

• 

• 

• • 

4 


In  this  Table  the  first  digit  of  oach  subscript  specifies  the 
position  of  the  source  oloment,  and  has  the  range  1 to  4,  soo  Fig.  1.  The 
second  digit  specifies  the  row  of  the  othor  element,  and  has  the  range 
1 to  7.  From  considerations  of  symmetry,  it  is  found  that  for  the  four 
source  positions  and  any  one  column  of  elements,  there  exist  16  different 

values  of  Z.  ..  Each  dot  in  Table  2 represents  a value  of  Z . , equal  to  one 

„ ‘J  . >J 

of  the  values  written  in,  as  givon  in  Tablo  3. 

TABLE  3 


Z33 

Z55 

Z17 

ZZ 

Z71 

Z22 

= 

Z66 

Z35 

= 

Z53 

Z11 

n 

Z77 

Z26 

- 

Z62 

Z12 

Z21 

= 

Z67 

= 

Z76 

Z13 

= 

Z31 

— 

Z57 

= 

Z75 

Z14 

Z41 

s 

Z47 

= 

Z7  4 

Z15 

Z51 

= 

Z37 

= 

Z73 

Z16 

= 

Z61 

- 

Z27 

= 

Z72 

Z23 

= 

Z32 

= 

Z56 

= 

Z65 

Z24 

= 

Z42 

= 

Z46 

= 

Z64 

Z25 

= 

Z52 

= 

Z36 

= 

Z63 

Z34 

= 

Z43 

= 

Z54 

= 

Z45 

There  are  68  columns  of  olomonts  in  the  transducer,  of  which 

34  v/ i 1 1 yiold  difforent  values  of  Z.  .,  by  symmotry.  The  total  number  of 

J 

different  values  of  Z.j  existing  for  the  trmsducor  is  therefore  16x34=544, 
In  the  general  case  of  a transducer  of  this  shape  and  symmotry 
having  N rows  and  M columns,  (M,N>2)  the  number  is  the  product  of 
(1+3+5+.. .N)  or  (2+4+6+. ..N)  for  N odd  or  oven  respectively,  and  (f4— 1 )/2 
or  f/2  for  M odd  or  oven,  respectively.  It  may  bo  noted  that  M and  M are 
not  interchangeable  in  thoso  expressions;  that  is  bocauso  the  vertical 


l 
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symmetry  (mirror)  and  the  horizontal  symmetry  (circular)  are  different. 

For  the  transducer  shown  in  Fie.  1,  in  addition  to  the  544 
different  values  of  Z.  . there  are  also  four  different  values  of  self- 
impedance Z.,  one  for  each  of  the  four  distinct  source  positions. 


10.  Computation  Req'ijrad 

It  is  desired  to  calculate  values  of  / p and  Z.  . for  the 
various  transducer  elements.  This  requires  bothm*  calculation  of  the 
spheroidal  functions,  which  are  not  tabulated  in  the  range  required,  as  well 
as  the  numerical  evaluation  of  certain  integrals.  A discussion  of  the  scopo 
of  this  work  is  given  in  the  following  sections. 


10a.  Angular  Functions 

In  the  present  case  the  transducer  is  11. 4\  in  diameter,  and  the 
correspond ing  frequency  parameter  c f or  the  spheroid  has  the  relatively 
large  value  of  32.8.  None  of  the  existing  tables  of  spheroidal  functions^ 
extend  to  such  high  values  of  c,  and  thus  the  required  angular  and  radial 
functions  must  be  computed. 

Flammer  gives  an  asymptotic  expansion  for  the  angular  functions 
that  is  valid  when  c|  is  large,  as  in  the  present  case  where  c-  = 19.68; 
see  Ref.  10,  pp  62-7,  equations  (3.2.9)  to  (8.2.42).  The  expansion  is  in 
terms  of  the  associated  Laguerre  polynomials,  and  should  give  an  easier 
computing  task  than  the  usual,  non-asymptot ic  expansion. 

We  require  S(nn(-32.8i,»p  for  going  from  0 to -0.65;  for  the 

lower  modes,  only  a few  values  of  will  be  required  to  define  the  function 

accurately  over  the  range;  for  the  higher  modes  more  values  of  will  be 

needed.  The  resulting  values  of  S (>1)  will  be  used  first  in  6.12  to 

0 mn  v 

evaluate  the  expansion  coefficients  specifying  the  boundary  condition  and 
then  again  in  7.3  to  obtain  the  pressure  at  a field  point. 

It  is  not  known  at  present  how  many  modes  will  be  required  to 
fit  to  the  boundary  conditions  with  reasonable  accuracy;  however,  since 
the  transducer  element  has  an  angular  width  of  some  5.4°  which  is  - 1/33 
of  the  180°  range  from  pole  to  pole  of  the  spheroid,  it  seems  probable  by 
considering  the  corresponding  Fourier  expansion  that  not  less  than  30 
harmonics  would  be  needed  to  provide  accuracy  of  the  order  of  10  °/o  . 

11  J.  A.  Stratton,  P.  M.  ’torse,  L.  J.  Fhu  and  R.  A.  Hutner,  Cll  iptic 
Cylinder  and  Spheroidal  ■’  ove  Functions.  John  iley  and  Sons  (1941). 

12  J.  Stratton,  P.  M.  Morse.  L.  J.  Chu,  J.  D.  C.  Little,  F.  J. 

Corhitb,  Suharoidal  ■■  ivo  Functions.  John  - iley  and  Sons  (1996). 
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The  norm  N of  the  annular  functions  is  required  in  6.10.  On 
mn  1 

p.  67  of  Ref.  10  an  asvmptotic  exoression  R for  this  quantity  is  discussed. 

' tnn  ' 


10b.  Radial  Functions 

The  pressure  is  needed  only  on  the  surface  of  the  spheroid,  thus 
only  one  vilue  of  the  radial  function,  and  one  of  its  derivative  need  be 
calculated  for  each  mode.  Since  £o=0.6,  and  c=32.8,  cJ0=19.68,  and  the 
latter  value  is  largo  enough  for  the  following  asymptotic  expression^  to 
be  used: 

if  exp  4 [c$"  2 (nrf1)  (10*1) 

The  real  part  of  this  radial  function  corresponds  to  the  pressure  component 
in  phase  with  the  velocity,  and  the  imaginary  part  to  that  out  of  phase. 

The  asymptotic  expression  for  the  derivative  of  the  radial 
function  can  be  obtained  directly  from  9.1  by  differentiating: 

^mn^  exp  i [c£-4f(n+1)  tt]  (10.2) 

ci 

Thus  from  10.1  and  10.2  the  ratio  of  the  radial  function  to  its 
derivative,  as  roquireJ  in  8.3  is 

r(3) 

— ^yr  = j/(ic.f-  1)  (10.3) 

R 

mn 


10c.  Expansion  Coefficients 


In  6.10  the  expansion  coefficients  a • for  a mode  mn  and  position 

mn t 

i of  the  transducer  element  are  given  in  terms  of  an  integral  defined 

in  6.12.  may  be  calculated  as  follows: 

1)  From  6.4  hj(»j)  is  calculated  over  the  range*]  =0  tov|=)/^,  whore  vj^ 
corresponds  to  the  half-hoipht  of  the  transducer.  The  same  function  h^(>\) 
will  be  used  for  all  modes. 

2)  S (h)  is  obtained  as  discussed  in  10a. 

mn  1 

3)  From  (1)  and  (2)  the  product  sjv  is  thon  formod  for  each 

mode  mn. 


4)  The  integral  in  6.12  is  required,  and  this  is  obtained  numerically 
for  each  mode  mn  over  the  ranees'*]  jg*  * = ^ to  4j  this  yields  4 mn  valuos. 

5)  Each  of  the  values  obtained  in  (4' is  multiplied  by  ^ sin  m 0^, 

13  Flammor  in  Ref.  10  gives  10.1  incorrectly  on  p.  67,  Eq.  (8.2.4R)  as  the 
limit  as2)-»c<^  but  correctly  on  p.  32,  Eq.  (4.1.17)  as  limit  when  cj  —too. 
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yieldino  4mn  values  of  L ..  The  factor  u can  hn  carried  through  as  a 
' s onn 1 

s^nbolj  it  will  go  out  when  the  expression  for  (8.2)  is  calculated. 

Having  found  those  values  of  1^.,  arrB-,i  *s  obtained  from  6.10, 

as  in  Section  10a. 


/ _V  | III*  I • I'M  I * 

R' ^ ) havinp  been  obtained  as  in  Section  10b,  and  N 

mn  o 

There  are  4mn  values  of  a_  . . 

mm 


mn 


1 0d . Calculation  of  Pressure  .it  Transducer  jo, 

Following  the  scheme  of  calculation  given  in  10a,  b and  c, 

21  P c ,0)  as  pivon  in  7.3  is  now  calculable  as  a function  of  angular 

nn  mn i 1 

position  (T),0)  on  the  spheroid,  using  as  many  modes  as  are  needed.  Z.j  is 
then  calculated  from  8.1,  evaluating  the  integral  numerically. 
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List  of  Symbols 


a . = numerical  coofficiont  for  mn  mode,  in  expansion  of  souno  field 

mm  ^h 

produced  by  driving  tho  i transducer  el  orient 

c = frequency  parameter,  = kd/2 

c = sound  velocity 

o 

d = intor  focal  distance  of  spheroid 

hjC^)  = coordinate  scale  factor 
i,j  = transducer  elements 

1^.  = integral  defined  in  (6.6) 

k = wave  number,  = (q/cQ 

■^mni  = nurTier*ca^  coefficient  in  expansion  of  boundary  condition 

m,n  = numbers  specifying  tho  vibrational  mode 

M,N  = integers  > 2 

N = norm  of  S ('l) 

mn  mn 

p.  = pressure  field  radiated  by  i element 

o . = pressure  of  mn  mode  radiated  by  i ^ element 

Kmni  r 

R^(-ic,  il)  = radial  function  of  the  3rd  kind,  of  order  m and  degree  n 
mn 

R^'  (-ic,i'|)=  derivative  of  R^(-ic,  if)  with  respect  to 
mn  mn 

S M = anpular  function  of  tho  1 kind,  of  order  m and  degree  n 
mn  1 13 

t = t imo 

u = velocity  amolitude  of  driven  transducer  element 

Uj»Uj  = velocity  of  i^  and  elements 
v^  = velocity  normal  to  spheroid  surface 
v = given  velocity  normal  to  spheroid  surface 

x,y,z  = rectangular  coordinates 

l.  . = mutual  impedance  between  elements  i and  j 

^ th 

2 = radiation  impedance  of  j element 

J 
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1 

^«1’^i2= 

nn  = 

X 

* \ = 
yio 

p 

0 

0i1^i2  " 

-Wi  = 

Y ■ 

Yi  * 


angular  coordinate 

angular  coordinates  defining  the  upper  and  lower  sides  of 
the  element 

angular  coordinate  dafinino  upper  side  of  element  in  position  1, 

Fig.  1 

wavelength 
radial  coordinate 
fluid  density 
angular  coordinate 

angular  coordinates  defining  the  vortical  sides  of  the  i^  element 

angular  coordinates  defining  the  vertical  sides  of  the 
i*^  elemont  when  i=1  to  4,  Fig.  1 

velocity  potential 

velocity  potential  caused  by  radiation  of  i^*1  element 
2 7T  x f requoncy 


ft 
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